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Abstract 

A geometrical model which captures the main ingredients governing 
atom-diatom collinear chemical reactions is proposed. This model is 
neither near-integrable nor hyperbolic, yet it is amenable to analysis 
using a combination of the recently developed tools for studying systems 
with steep potentials and the study of the phase space structure near 
a center-saddle equilibrium. The nontrivial dependence of the reaction 
rates on parameters, initial conditions and energy is thus qualitatively 
explained. Conditions under which the phase space transition state 
theory assumptions are satisfied and conditions under which these fail 
are derived. 

1 Introduction 

The study of classical, semi-classical and quantum chemical reactions on a 
molecular level has a rich history [UEIEIIHEIEIEIIEJE]- ^ n these models, the 
full Hamiltonian is averaged over the fast motion of the electrons, where each 
electron is assumed to be fixed at a specific quantum energy level (the adia- 
batic approximation) [U [2j |3]. Such computations produce effective potential 
energy surfaces (PES) that govern the slow motion of the nuclei. The resulting 
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Hamiltonians correspond to the "Born-Oppenheimer" approximations. Quasi- 
classical computation^] that employ these Hamiltonians provide surprisingly 
good approximations to the quantum calculations, hence classical models are 
extensively studied by chemists, see pQ [21 [3J and references therein. 

Dishearteningly, the resulting nuclei motion of even the most basic, classical 
"Elementary" (bimolecular) reactions, that are "at the heart of chemistry" [3J, 
is not well understood. Numerical simulations of the nuclei dynamics exhibit 
sensitive dependence of the trajectories on initial conditions and parameters. 
Moreover, one finds that macroscopic observables, like reaction rates, have an 
intricate dependence on the parameters and energy, see e.g. [H 121 El El IB] • 

In contrast, the popular transition state theory provides an appealing intu- 
itive view of this motion and leads to explicit formulae relating the microscopic 
nuclei kinetics to the macroscopic reaction rates. This theory assumes that the 
PES is separable to a one dimensional potential ("along the reaction path") 
and a potential well in all the other modes of motion (the bath of oscillators). 
Hence, the reaction according to this theory is described by a product of a 
one degree of freedom (thus integrable) system and oscillators. This appeal- 
ing phenomenological model is too simple: it does not describe the observed 
complex dependence of the nuclei motion on initial conditions and energy 

E 121 12 El mi. 

The next level of approximation, by which the non oscillatory part of the 
potential is two-dimensional, is the subject of our paper. Here, the transla- 
tional and vibrational energies of the atom and diatom are coupled to each 
other yet are separable from all the other modes of motion (e.g. of those 
associated with the bending and the rotational energies). Such a decoupling 
occurs, for example, when the initial configuration is collinear and the angular 
momentum is zero [U E]. This separability assumption is widely used in the- 
oretical investigations of chemical reactions [HEIEIIIIEIEIEIIHIE]- Indeed, 
the first examinations of the principles underlying the transition state theory 
were conducted by investigating such models [U E] • 

Employing such a restrictive reduction)^] still leaves us with a two degree 
of freedom system which may admit chaotic behavior (in contrast with the 
transition state theory reduction which leaves us with integrable dynamics). 



1 In these computations the atomic motion is found from the classical dynamics dictated 
by the PES. The quantization enters through choosing initial quantised ensembles (usually 
only in vibrational and rotational energies) and through binning of the product states to 
quantized energies. 

2 By which the bath of oscillators is considered to be separable from the reaction dynamics. 
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One source for such complicated behavior is associated with the existence of a 
reaction barrier - a saddle point of the PE^J The existence of such a barrier 
is the main ingredient needed for relating the classical transition state theory 
to the analysis of the phase space structure of two degrees of freedom systems 
[U H21 EES HH HH1 HE] ■ Normal form analysis of the local phase-space structure 
near the PES bottlenecks (local unstable extremal points of the PES) relates 
these two approaches and leads to accurate calculations of the minimal flux 
through them [T2], \TU\ [11]. The local picture near the bottlenecks does not 
reveal the complexity of the motion. This complexity is revealed only when 
the global structure of the reaction tubes - the phase space regions that pass 
through the bottlenecks - is calculated |S1 El H21 H31 [HI HH UHl HHJ- 

Previous works on the global features of these tubes have either employed 
near-integrable techniques [HJ [201 HU [21] or numerical integrations Q2J 
IT3| HU [T71 HH [19] . In these works it is demonstrated that the structure of the 
stable and unstable manifolds of unstable periodic orbits (the Periodic Orbits 
Dividing Surfaces PODS [5j or, similarly, the stable and unstable manifolds of 
the Normally Hyperbolic Invariant Manifolds NHIM [12[[T8]) determines the 
reaction rates. In particular, when these asymptotic surfaces intersect each 
other the reaction rates depend on both the intersection pattern and the local 
flux near the saddles. In such cases the predictions of the transition state 
theory (and its local minimal flux variants) must be modified. These insights 
are employed to study numerically the asymptotic surfaces to NHIMs of high 
dimensional problems p2J [17] . 

Here we introduce a new geometrical model for the two dimensional PES. 
The PES in the reaction region is a approximated by a sum of a quadratic 
form with a saddle point and a smooth potential which is close to zero within 
a corner region and increases sharply at the corners' boundaries (see Figs[l]3j2]). 
We analyze this model by combining the theory of smooth Hamiltonians with 
impacts j22j |23l EH [25], the recent generalization of [261 EZ] to the impact 
case [28J, and the theory of homoclinic loops to a saddle-center [291 USB I3T] . 
The analysis provides qualitative understanding regarding the global structure 
of the dividing surfaces in reactions with one rank-1 saddle point. Indeed, 
we identify a basic mechanism which explains why and when the dividing 
surfaces of the periodic orbits that have energies slightly above the barrier 

3 Some reactions have potentials that also admit stable triatomic states (indirect reac- 
tions), many reactions have a single steady unstable triatomic configuration (direct reac- 
tions), and in some exceptional cases there are no steady triatomic configurations at all 
[3] . Most recently, models with rank-k saddles have been considered as well [lOl [11] . 
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energy are especially complicated or especially simple (namely do not intersect 
each other). Moreover, we provide a mechanism for the emergence of stable 
triatomic cyclic motion. The applicability of these geometrical insights to more 
accurate models of the PES and to higher dimensional settings (e.g. when the 
two degrees of freedom dynamics are weakly coupled to a bath of oscillators) 
is under current study. 

The paper is ordered as follows. In section 2 we recall that mass-scaled 
Jacobi coordinates bring the collinear triatomic reaction's Hamiltonian to a 
two-degree of freedom Hamiltonian in the standard mechanical form. In sec- 
tion 3 we introduce the geometrical potential function and in section 4 we 
provide some basic observations regarding its properties and its relation to 
other potentials. In sections [5] we analyze the dynamics in this model under 
specific geometrical conditions, proving that both chaotic and stable periodic 
triatomic motions emerge. In section 6 we show that complicated dynamics 
appear at some parameter ranges and simple dynamics at others. In section 7 
we discuss the relation of these results to transition state theory and propose 
some possible extensions. 



2 Collinear atom-diatom reactions 

Some of the geometrical characteristics of the PES describing collinear tri- 
atomic reactions may be inferred from general considerations that are common 
to all such reactions. These characteristics, as described next, motivate our 
construction of the simplified geometrical model for the PES. 

Consider the triatomic reaction A+BC — > AB+C. Namely, here we always 
consider the atom-diatom case. We denote this reaction by the standard short- 
hand notation AB + C . An effective Hamiltonian for the molecular interaction 
is of the form 

H(r, P )= £ ~ + V{r) 

i£{A,B,C} 1 

where r = (r^, r^, r^) denotes the positions of the atoms A,B,C in a given 
inertial frame and {Mi} i( -^ A B c y denote the masses of the atoms. Under stan- 
dard conditions V depends only on the relative positions of these atoms. This 
9 d.o.f. system simplifies when collinearity is assumed [HIS]. This assumption 
implies that the relative positions may be expressed in terms of two scalars 
r i = ( r A — t b) ■ e, r 2 = {r B — r c ) • e where e is a unit vector aligned with the 
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molecules, namely V(r) = V(ri,r 2 ). Moreover, since at small distances the 
atoms are strongly repelling, V becomes large along the rays = 0, see Fig. 
[T^i. The kinetic energy term in these new coordinates is non-diagonal and has 
a mass dependent quadratic form. A mass-weighted coordinate system (the 
Jacobi coordinates, see [2] for formulation and references) brings the system 
to the standard mechanical Hamiltonian form of a unit mass particle moving 
in the potential field V r (qi, q?): 

2 2 

H(qi, P i) = ^ + ^ + V r ( qi ,q 2 ). (2.1) 

Here, ((71,(72) are the reaction coordinates: 

qi(ri,r 2 ) = a r 1 + b r 2 cos (3, q 2 (r 2 ) = b r 2 sin (3, (2.2) 
the scaling coefficients a, b depend only on the mass of the atoms: 



M A (M B + M c ) t / M C (M B + Ma) 



Ma + Mb + Mc' y Ma + Mb + Mc 1 

and the mass dependent "skew-angle" (3 is defined by: 



(2.3) 



M A M C ( . 

p = arccos \ / — — — r— — 7~rw- (2.4) 

For the H 2 + H reaction, (3 = 60° (see FigjlJ)), for heavy-light-heavy interac- 
tions (3 is small (e.g. for IH + / one finds f3 = 7°) whereas light-heavy-light 
interactions lead to j5 ~ 90°, see [TJ EJ El [2] for the corresponding figures and 
references therein. 

The short range repulsion of the atoms implies that the motion in the 
configuration space is confined by the rays r\ = 0, r 2 > (namely q 2 = 
gitan/3,g 2 > 0), and r 2 = 0, ri > (namely q 2 — 0, qi > 0). This region, 
a two-dimensional wedge in the (qi, q 2 ) plane, is called hereafter the corner re- 
gion or the /3-wedge, see Fig[T]D. The dynamics within the /3-wedge depend on 
the particular form of the potential V r (qi,q 2 ). We will be mostly considering 
potentials that have a single extremal point which is a saddle. 

In the scaled coordinates, chemical reactions are represented by trajecto- 



ries of the Hamiltonian (2.1). Reactant states (BC molecules and far away 
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A atoms) correspond to configurations with a large qi and bounded q 2 . We 
thus say that such configurations belong to the reactants channel. Similarly, 
product states (AB molecules and far away C atoms) correspond to configu- 
rations in the products channel with bounded (gitan/3 — q 2 ) and large q±,q 2 . 
In the reactants (products) channel the potential is well approximated by a 
one dimensional BC (AB) diatomic potential. The corner region, where both 
qi and q 2 are bounded, and where all the potential saddle points are located, 
is called the reaction region. In symmetric cases with a single saddle point 
the saddle point is located on the bisector - the potential symmetry line. In 
non-symmetric cases, the barrier location is called 'early' (respectively 'late') 
if it is closer to the reactants (respectively products) channel. 

3 The geometrical potential function 

The above description of the geometrical properties of the potential V r (qi, q 2 ) 
is independent of the details of the reaction. Below, we propose a specific form 
for a potential V r (qi,q 2 ) with parameters that have a transparent geometrical 
meaning. This potential may be viewed as a local geometrical approximation 
in the interaction zone to any other potential surface. The advantage of using 
this new formulation becomes apparent - it allows for rigorous analysis of the 
model and for qualitative understanding of the dynamics. 

Consider Hamiltonians with geometrical potentials of the form: 

p 2 

H(q,p; a,b,c,e) = — + aV a (q) + bV b (q,e) + cV f 

ar field 

(q). (3.1) 

Here H\, = H(q,p; b, a = c = 0) is a billiard-like system limiting to a billiard 
in a /3-wedge (see more details below). H int = H(q,p;a, b = c = 0) is an 
integrable system where the potential V a (q) has a single saddle point in the 
corner region (we will soon fix V a (q) to be a quadratic potential and remark 
on the effects of higher order terms when applicable). The far-field potential 
cVf ar fieid{q) and its derivatives are small in the reaction region (yet are large 
away from the corner region). All together, aV a (q) corresponds to the normal 
form of the potential near the saddle point, bVb(q,e) corresponds to the di- 
atomic repulsion terms and cVf ar f ie id(q) handles the remainder terms near the 
saddle and the reactants and products channels away from the corner region. 

By billiard-like system we mean that the potential Vb(q, e) is a steep po- 
tential [32]: the level set of this potential at, say, = 1/2, limits, as e — > 
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to some billiard-like domain and Vb(q,e) — > for all q in the interior of this 



domain (see [27] for a precise definition and Eqs. (3.4), (3. 5) for examples). 
Here, the /3-wedge is the billiard domain. For example, in Fig Vp, the red level 
curves are identified with the level curves of Vf,(q,e). For small a values (or 
equivalently, at energy levels that are much larger than a), the smooth part 
of the potential at the corner region may be neglected and the motion in the 
/3 -wedge is billiard-like [26|l27]. 

For non-negligibl^ja, for sufficiently small (e, c), at a fixed positive energy 
level, the system H e = H(q,p; a, b, c, e) has trajectories that closely follow the 
smooth integrable dynamics H int until they reach the walls defined by H(<?, £), 
reflect according to the billiard law, and continue with the smooth dynamics. 
Indeed, in the limit e — > 0, the motion is described by a Hamiltonian system 
with impacts [221 EH EU ES] . Here, we study the dynamics of such impact sys- 
tems when their smooth integrable part (Hi nt ) has invariant hyperbolic subsets 
(here, Lyapunov periodic orbits near a saddle-center [331 [2S1 [30]). We prove 
that for a range of parameter values the reflections of the stable and unsta- 
ble manifolds of this set from the billiards' boundary give rise to complicated 



homoclinic behavior (Theorem 5.2) and to stable, recurrent triatomic states 



(Theorem 5.5). We also show that at a different range of parameter values 



the manifolds reflect to infinity via the reactants and products channels and 



no recurrent motion near the saddle is possible (Theorem 6.1 and Conjecture 



6.5). For this latter range of parameters transition state theory is expected to 
be valid. 

Finally, the results established for the limiting impact flow at e = c = are 
shown to be valid at small e, c. The persistence for small c values follows, as 
usual, from the robust character of these results under smooth perturbations. 
The persistence for small e follows from the recent extension of [27J to smooth 
Hamiltonians that limit to impact systems [28]. 

3.1 The simplest form of the geometrical model 

The assumption that the potential has a single saddle point in the corner region 



implies that the potential aV a of (3.1) is of the form 

aV a (q) = \{q~ q s ) T A(q - q s ) + 0((q - q s ) 3 ) (3.2) 



*Yet bounded, so that along the corner boundaries the diatomic repulsion dominates. 
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where A is a symmetric 2x2 matrix (A T = A ) with eigenvalues (u 2 , — A 2 ). 
Additionally, since the potential saddle point q s = (g ljS , q 2tS ) is assumed to be 
located within the corner region, it satisfies gi )S > 0, q2. s < <?i.stan/3. Conse- 
quently (see section [1]), the Hamiltonian flow has a saddle-center equilibrium 
at P= (q,p) = (q s ,0). 



Replacing aV a (q) in Eq. (3.1) by its quadratic approximation we obtain 



the simplest form of the geometrical model: 



JT 1 

H(q,p;b,c,e) = — + -(q - q s ) T A(q - q s ) + bV b (q;e) + cV farfidd (q). (3.3) 

Recall that the bV b term corresponds to a billiard-like potential in the corner, 
namely it satisfies conditions I- IV of [27J. We may, for example, consider 
here either an exponential potential (as in the repulsion associated with a 
Morse potential) or a power-law potential (as in the Pauli repulsion term of a 
Lennard- Jones interaction) : 

bV b (q,e) =bV exp (q,e) = bexp(-q 2 /e) + 6exp((g 2 - <?i tan/3)/e) (3.4) 



or 



%<)=«UH,'W = jt ^ , kez+. (3.5) 



We study Eq. (3.3) in the small c regime, namely, we assume that the 
far-field term is small. We expect this approximation to be valid only in the 
reaction region. We thus study the dynamics in a bounded corner region of 
the configuration space. The simplest model provides adequate approximation 



to (3.1) in this bounded region if the far-field potential and all its derivatives 
are small there and, additionally, aV a (q) is well approximated by its quadratic 
approximation. This latter part of the assumption may be relaxed in the future 
by including higher order terms of the normal form near the saddle-center point 
of the Hamiltonian. 

The motion for a, b ^ is found by analyzing first the singular impact 
limit (i.e. when e — > +0), hereafter called the limit system. Then, the recent 
persistence results of [28] (generalizing [261 EZ]) an d standard perturbation 
theory allow us to establish that similar behavior persists for sufficiently small 
e and c. 
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4 The phase-space structure of the limit sys- 
tem and its dependence on parameters 

In the limit e — > +0 the behavior of ( |3.3[ ) in the interior of the corner region 
is governed by the linear integrable system corresponding to the quadratic 
Hamiltonian. The limit system is defined as this smooth integrable motion 
inside the corner together with reflections from the corners' boundaries. 



In section 4.1 we recall the integrable phase space structure of the quadratic 



Hamiltonian in the saddle-center case. In section 14.21 we recall the reflection 



law from the corners' boundaries. In section 4.3 we explain how the parameters 



of the simplest geometrical model (3.3) may be extracted from a general PES 



and motivate our choice of particular ranges of geometrical parameter values 
in sections I 



4.1 The integrable structure near the saddle-center 

Below, we first rotate the (q, p) coordinate system so that the quadratic Hamil- 



tonian in (3.3) becomes separable (more generally, one brings the integrable 
part of (3.1) to its normal form near the saddle-center point). We then de- 
fine the corner region in the rotated coordinate system (w,t>), see Fig [2] We 
show that in the rotated system the projections of the stable and unstable 
manifolds of the saddle-center fixed point and of the nearby periodic orbits are 
easily found. We end this subsection by defining two constants of motion and 
identifying reacting and non reacting trajectories near the saddle-center point 
in terms of the values of these constants of motion. 

Denote the normalized eigenvector of A which corresponds to u 2 (respec- 
tively (—A 2 )) by U\ (respectively 1/2). Recall that U1-LU2. A small computa- 
tion shows that the eigenvalues of the saddle-center point of the Hamiltonian 
flow, P = (q s ,0), are ±iu,±X. The Hamiltonian flow has a two-dimensional 
real invariant subspace E c (center plane) corresponding to the eigenvalues 
±iu : E c = {(q,p)\(q,p) = (rUi,sUi), r,s G E}. The center plane projection 
onto the configuration space is one dimensional, along the eigenvectoij^ U\. 
Similarly, the Hamiltonian stable and unstable subspaces corresponding to the 
real eigenvalues ±A are expressed in terms of the matrix eigenvector U 2 : E u = 
{(q,p)\&P) = r(U 2 ,XU 2 ), r e R},E" = {(q,p)\&P) = r(U 2 ,-XU 2 ), r e R}. 



If higher order terms of the normal form are included in aV a , these statements apply to 
the tangent planes of the corresponding manifolds. 
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Thus, the projection of both the stable and unstable subspaces onto the config- 
uration space is along 11% - the eigenvector of A corresponding to the negative 
eigenvalue (—A 2 ). The eigenvectors directions should not be confusecj^jwith the 
zero level-lines of V a . Since the potential V a has a saddle point at q = q s , there 
are two directions, Q ± = [/i + a^L^, along which V a vanishes V a (q s + tQ ± ) = 
for all t, and are finite and non- vanishing (see Fig [2]). 

To simplify further calculations it is convenient to transform the Hamilto- 
nian to its linear normal form. Let us denote by (cos#, — s\.n9) T the unitary 
eigenvector Ui, assuming, with no loss of generality, that — n/2 < 9 < ir/2. 
Then (sin^, cos6*) T is the unitary eigenvector U 2 - A useful standard observa- 
tion for natural mechanical Hamiltonians is that rotations of the configuration 
space can be easily incorporated into the Hamiltonian. Indeed, defining the 
standard rotation matrix: 

Re=( C °^l ~™$) (4.1) 



sin(6 l ) cos(#) 

and making the symplectic transformation: (q,p) —> (u,v) = (Rgq, R$p), the 



Hamiltonian (3.3) becomes: 

v 2 1, „ , 

(4.2) 



H(u, v; b, c, e) = —+-(u-R e q s ) T R AR_ e (u-R e q s )+bV b (R„ e u, e)+cV farfie i d (R- 



where the quadratic part is diagonal: 

-A 



R e AR_ e =[ U " ° ) . (4.3) 



The integrable part of the Hamiltonian (3.3) becomes (where u s = Req s ) 



H lm = + y K - u,,) 2 - ^(u 2 - u 2 , s ) 2 . (4.4) 

Recall that the quadratic (or more generally the integrable) approximation 
inside the corner is expected to hold only in the reaction zone where the far- 
field contribution is small. We thus define a bounded corner region in the 
(mi,m 2 ) plane, Al, and study the dynamics in this region (see Fig [2]). The 



6 Thcse simple observations are the source of much confusion as we wrongly tend to 
confuse the level sets of V a in the configuration space q with the phase space plots in the 



space spanned by the stable and unstable directions (see Figs. 2p ) 
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lower and upper boundaries of Al are the two rays that emanate from the 
origin and are aligned with the vectors (cos 9, sin 9) (lower boundary), and 
(cos(/3 + 0),sin(/3 + 9)) (upper boundary). This (3— wedge is then intersected 
by the square [-L, L] x [-L, L], where L > max(ui )S , U2 )S )- Since (3 G (0, 7r/2), 
we have (3 + 9 G (— 7r/2, 7r) and two cases appear. For (3 + 9 < n/2 (see Fig. 

A L = {«i G [0, L],tt 2 G [«i tan 0, min-{/ui tan(/3 + 9),L}]} (4.5) 

whereas for (3 + 9 > n/2 (see Fig. [2|d): 

Al = {ut G [0,L],« 2 G [«itan0,L]} U {u x G [cotan(/3 + 9),L],u 2 G [0,L]}. 

(4.6) 

Finally, we describe in geometrical terms the well known structure of the 
linear flow in Al- The motion under this linear flow occurs on surfaces defined 
by the joint levels of Hn n and the action Ix- 

2 2 

Ix(ux, ui) = y + y(«i - u hs f. (4.7) 

The action I\ is the constant of motion associated with the oscillatory motion. 
The other constant of motion: 

v 2 X 2 

D 2 (u 2 , V2 ) = -±- —{u 2 - u 2yS ) 2 = H lin - h (4.8) 

determines the hyperbolic motion in the (u 2 , t>2)-plane. The surface on which 
the motion occurs, {(u,v)\Hn n (u,v) = h,Ix(u,v) = k}, is composed, for k > 
and h ^ k, of two disconnected 2-dimensional cylinders: the direct product of 
an ellipse in the (ux, f i) plane and two branches of a hyperbola in the (u 2 , v 2 ) 
plane. The sign of D 2 = h — k determines the nature of the hyperbolic motion. 
For negative D 2 the hyperbola branches are directed sideways, one having 
positive u 2 — u 2)S and the other having negative u 2 — u 2 ,s> see the shaded 
region in Fig. [3} Trajectories belonging to level sets with D 2 < do not 
cross the 3-plane u 2 = u 2tS : they approach it and then return to the same 
side. Keeping in mind the chemical origin of our model, we shall say that the 
motion corresponding to such trajectories occurs "without reaction". On the 
other hand, if D 2 is positive, branches of the hyperbola extend horizontally 
along the full u 2 axis. All the trajectories that belong to level sets with D 2 > 
cross the surface u 2 = u 2yS . The upper branch of the hyperbola corresponds to 
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orbits with a monotonically increasing u 2 ("reactants to products"), whereas 
the lower one to monotonically decreasing u 2 ("products to reactants"). We 
thus say that such trajectories "realize the reaction". 

The level sets on which Hn n — h — k — I\ (so D 2 — 0), are the singular level 
sets. These sets separate the two types of motion (with vs. without reaction). 
Each such singular level set contains a normally hyperbolic Lyapunov periodic 
orbit 7^ = {Ii = h > 0,U2 = u 2jS ,v 2 = 0} belonging to the center plane E° 
along with its stable and unstable manifolds iy s,M (7/ l )(each being a straight 
cylinder). At h = this singular level set contains only the saddle-center point 
P and its stable and unstable manifolds W S ' U (P = 7 )(each being a straight- 
line). The projection of these local stable and unstable manifolds (W S,U (P)) 
onto the configuration space is a straight line, aligned with the vector U 2 . This 
projected line is divided by the saddle point u s into two rays: the extension 
of the lower ray, the projection of Wl ,u (P), intersects the lower boundary of 
the corner, whereas the extension of the upper ray, the projection of W+ U (P), 
intersects (if tan(# + /3) > 0) the upper boundary of the corner (see right panels 
of Fig. [2]). The projection onto the configuration space of the stable and 
unstable manifolds of a Lyapunov orbit 7^ with h > appears as a collection 
of many oscillatory orbits that are centered around this line. The projection 
of the Lyapunov periodic orbit 7^ lies within the corner region provided h is 
smaller than /i crit _ 7 : 



h CTit _^ = min <^ — (u 2)S / tan(/3 + 9) - u 1>s ) 2 , — (u 2;S / tan0 - u ltS ) 2 \ . (4.9) 



For large h values the impacts destroy these periodic orbits. Hereafter we 
always consider energies that are close to the saddle point energy and are thus 
strictly smaller than h CTit _^. 

Summarizing, for any given h, the energy surface Hu n = h is foliated by the 
levels of I\ into the level sets on which the motion occurs. For h < 0, all these 
level sets have negative D 2 , hence, these do not cross the surface u 2 = u 2 , s 
and no reaction may occur there. When h > the energy surface contains, 
additionally, a Lyapunov saddle periodic orbit with action I\{^h) = h and level 
sets with < h < h(lh)- These level sets have positive D 2 and thus orbits 
belonging to these correspond to reacting trajectories, namely, these orbits 
belong to the reaction tubes. 

See [3l] for a detailed explanation of the very similar analogous geometry 
of the energy surfaces near saddle-center-center-.. -center points in the higher 
dimensional settings and when higher order terms of the non-resonant normal 
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forms are incorporated. Here we concentrate on the two degrees-of-freedom 
case: the higher dimensional saddle-multi-center case with impacts may be 
studied similarly, leading to more complicated dynamics (involving, for exam- 
ple, homoclinic orbits to invariant tori as in |35j). 

4.2 The impacts 

In the rest of this paper we study how the standard integrable behavior changes 
when the trajectories of the linear system are reflected from the walls of the 
billiard corner. Here we recall the reflection law from the lower and upper 
boundaries of the corner region. 

The unit vector defining the upper boundary of the corner in the (u\,u 2 ) 
plane is sp+e = (cos(/3 + 9),sm((3 + 9)) T and its inward normal is n^+g = 
(sin(/3 + 9), — cos(/3 + 0)) T . The resulting reflection law is: 

y = (v h v 2 )^ , , 

((t>i cos 2 (/? + #) + t> 2 sin 2(/3 + 0),^ sin 2 + - t> 2 cos2(/3 + #))). y ' 

It is defined for velocities that exit the corner region, namely those satisfying 
(v, np+g) < 0. Similarly, for the lower ray, sg = (cos 9, sin 9) T and rig = 
(— sin 9, cos9) T , so the reflection law becomes: 

v=(vi,v 2 )H. , , 

(vi cos(2fl) + v 2 sin(20), v 1 sin(2fl) - v 2 cos{29)). y 1 



Remark 4.1 The reflection law preserves energy, namely the integral H . How- 
ever, reflections from the lower boundary (respectively upper boundary) do not 
preserve the integrals I\ and D 2 whenever 9 ^ 0,7r/2 (respectively whenever 
P + 9^0,ir/2). 

The change in the integrals by the reflections leads to the non-trivial behavior 
of the impact system. 

4.3 The geometrical parameters 

We show in sections |5]j6] that the dynamics of the limiting system depend in an 
essential way on the location and orientation of the saddle point with respect 
to the corner. The coordinates of the saddle point u s = Rgq s , the angle (3 of 
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the corner, the ratio u/X and the angle 9 between the eigenvector U\ and the 
gi-axis all matter in an essential way in determining the dynamics. 

These geometrical parameters, (/3, u s , L, 9, u, A, b, e), may be extracted nu- 
merically from any potential surface describing triatomic reactions. First, one 
finds the saddle point location q s , and linearizes the vector field at this point 
to obtain the matrix A. The angle 9 may be calculated via A's entries a^: 

tan# = J ~ 



y/(au -a 22 ) 2 + 4af 2 + a n - a 22 
The location of the saddle point in the normal coordinates is then found via 



u s = Req s (see 4.1). The estimates of b,e may be extracted from the diatomic 
potentials: these parameters are determined by the form of the strong atomic 
repulsion at short distances. Estimating the linear zone range L may be more 



delicate. In principle, comparing the approximation (4.4) to the numerical 
potential energy surface provides the range of validity of the linear approxi- 
mation. If it appears to be too small, it is possible to extend our theory by 
including higher order terms of the integrable normal form near the saddle 
center (as explained in detail, in another context, for the high dimensional 
chemical reaction settings in [M]). The application of this scheme to concrete 
reactions is under current study. 

A detailed classification of all possible dynamical behaviors for various 
(/3, u s , L, 9, uj, A, b, e) and h in the symmetric and asymmetric cases is beyond 
the scope of the current paper (moreover, it is hardly possible at all). In section 
5] we analyze the nearly perpendicular behavior (9 nearly zero or equivalently, 
9 + (3 close to zero). In section [6] we provide rough classification for the depen- 
dence of the manifolds' geometry on u/X. 



5 The nearly perpendicular dynamics 

Recall that the projections of W^ U (P), the lower branches of the local stable 
and unstable manifolds of the saddle-center point P, onto the configuration 
space are directed downward along the vector [7 2 . We show next that when 
U2 is close to being perpendicular to the lower boundary of the corner (so 9 
is nearly zero), near integrable behavior occurs. We establish first that when 
9 = e = c = the limit motion in some region containing the saddle-center 
point is integrable. We then prove that when 9 7^ the picture changes 
dramatically, leading sometimes to chaotic dynamics and sometimes to stable 
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triatomic periodic motion. The same results apply to the upper branches of 
the manifolds when (/3 + 9) is small. These two cases may arise, for example, 
in light-heavy-light reactions (/3 is close to 7r/2) with late (small 9) or early 
(small ((3 + 9)) barriers. 



5.1 Integrable behavior of the perpendicular limit sys- 
tem 



When 9 = e = c = and the energy h is smaller than /i cr jt_ 7 (see (4.9 )), homo- 
clinic loops are created by Wt ,u (-yh)- Indeed, as shown below, the projection 
of the stable and unstable manifolds of the center-saddle point P (respectively 
of the Lyapunov orbit 7/J is a straight line (a cylinder) which is perpendicular 
to the lower boundary of the billiard corner. Thus, after one reflection these 
manifolds coincide, see Fig. |4j 

Proposition 5.1 Consider the limit system (e — c — 0) at 9 — at an energy 
level h e [0, h crit - 1 ). The lower branch of the unstable manifold of 7^ coincides 
after reflection with the lower branch of its stable manifold, forming a family 
of homoclinic orbits to 7^. The flow of this limit system near the surface of 
homoclinic orbits is locally integrable: all nearby orbits belonging to the same 
energy surface h either belong to invariant tori or leave the homoclinic loop 
region after one round to the U2 > u 2 ^ s region. 

Proof. We first prove the existence of a homoclinic orbit for h = 0, where 
7^ =0 = P. For the case 9 = the projection of the stable /unstable manifolds 
of P onto the configuration space is a line perpendicular to the line u 2 = 0. 
Indeed, the linear system 

Ul = V!, U 2 = V 2 , Vl = -W 2 («i - Ux iS ), V 2 = \ 2 {U 2 - U2, s ) (5.1) 

has the flow 

V? 



Ui(t) — Ui iS = (u^ — Ui )S ) cos(wt) H — - sin(u;t) 



id 



v\(t) = —lu(ux — ui s ) sm(ut) + v\ cos(u;t), 

v° 

u 2 (t) - m 2 , s = («2 _ M 2,s) cosh(At) + -y sinh(At), 

A 

v 2 (t) = A(w!] — u 2tS ) sinh(At) + cosh(Ai). 



(5.2) 



So, the stable W S (P) and unstable W U (P) one-dimensional manifolds of the 
equilibrium P are the straight lines {u\ = ui tS , v\ = 0, u 2 = u 2tS — v 2 /\ } 



15 



and {u\ = Ui tS , V\ = 0, u 2 = w 2jS + respectively. Each straight line 

is divided by the point P into two rays W+ S (P) and W^ ,S (P). The lower 
rays Wl ,s {P) intersect the wall u 2 = (and the other rays intersect either the 
upper box boundary or the upper corner boundary). The stable and unstable 
rays that hit the lower wall intersect it at two different phase space points 
m s = (ui tS , 0, 0, Xu 2tS ) and m u = 0, 0, — Au 2jS ), respectively (recall that 
Ui jS ,U2, s axe both positive for sufficiently small 9). We choose, near each of 
these two points, sufficiently small 3-dimensional cross-sections to trajectories 
in the phase space: N s,u = {(u, v)\u2 = 0, \ \(u, v) — m S)U || < 5}. Each of these 
cross-sections is foliated into 2-disks N% = N s n{H = h}, N% = N U D{H = h}. 
As coordinates on the disk N% we take (u%,Vi), since the third coordinate t> 2 
on N s is expressed from H : v 2 = \/2h — (vf + u 2 (ui — Mi, s ) 2 ) + A 2 (u 2)S ) 2 = 
y/2h - 2/i + A 2 (« 2 , s ) 2 . The same coordinates (ui,vi) work for iV^, where the 
-0 2 -coordinate has the same form as v 2 but with a " — " sign in front of the 
root. If a trajectory of the linear flow hits N%, it is transformed to due to 



the reflection law (see (4.11)): the coordinates Mi,M 2 ,fi at the incidence point 
remain the same, but t> 2 = ?i 2 changes its sign. This means that the reflection 
law defines the symplectic global map (gluing map) Sh '■ N% — > N% as follows: 
U\ = iii, V\ = V\. In particular, m u is transformed to m s : we get a homoclinic 
orbit T to P. 

More generally, since 6 = 0, for any trajectory hittingTjthe 3-plane u 2 = 0, 
Sh simply changes the sign of v 2 . In particular, Sh preserves the energy and the 
integrals of motion l\ = Ii(ui,Vi) = h(ui,Vi) = k and -D 2 (w 2 ,f 2 ) = h — k = 
L>2- Thus, trajectories belonging to the level set h(u, v) = k, Hu n (u,v) = h 
that hit the lower boundary remain, after reflection, on the same level set with 
I\(u,v) — I\ = k and D 2 = h — k = D 2 : the reflection just changes their 
relative position along the same level set of the integrable linear system (see 
Fig§. 

Thus, for h G (0, /i CT .j i _ 7 ), the stable (I\ = h,D 2 = 0) and unstable {I\ = 
h, L>2 = 0) manifolds of the Lyapunov periodic orbit 7^ are cylinders that 
coincide after one reflection. Moreover, since (3 < 7r/2, these cylinders hit the 
3-plane u 2 = at u\ values that are bounded away from 0, namely bounded 
away from the corner. 

The dynamics for trajectories near the homoclinic loops are also fairly 
simple. In particular, trajectories lying on the given level set H = h, I\ = k 
are projected onto the plane (w 2 , v 2 ) into one of the two hyperbola branches 



7 i.e. trajectories arriving to this section with v% < 0. 
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of D 2 = h — k. If h — k < and u 2 < u 2 , s , then trajectories starting on iV 5 
move towards P and then return to the cross-section N u . In other words, the 
reflection law glues this hyperbola-like branch into a closed curve (see Fig. |4j 
the green triangle on N s is mapped by the flow to the green triangle on N u and 
then reflects to the black triangle, namely, back to N s ). Since such orbits also 
belong, in the («i,i>i) plane, to the closed curve I\ = k, these orbits belong, 
topologically, to an invariant torus in the phase space. Similarly, trajectories 
for which D 2 = h — k > belong to two disjoint cylinders (direct product of the 
circle I± = k in the (u\, v i)-plane and the two branches of the corresponding 
hyperbola in the (u 2 , t>2)-plane). Here, the gluing map defined by the reflection 
law glues the two end circles (those intersecting the section u 2 = 0) of these 
cylinders. Trajectories move along one of the cylinders towards the section 
u 2 = and after the first reflection escape along the second cylinder to the 
region u 2 > u 2iS . Their global dynamics, after they pass the saddle point, 
depend on ((3,u s ,co, X) as discussed in section [6j ■ 

We conclude that for 6 = e = c = the motion near T in the region u 2 < 
u 2]S is indeed simple. For negative energy h, the motion neai^jthe homoclinic 
loop occurs within the solid torus < I\ < II, D 2 = h — 1\ < 0. The motion 
is quasi-periodic if the rotation number on the related invariant 2-torus is 
irrational and it is periodic if this number is rational Inside this solid torus 
there is also a unique elliptic periodic orbit Ii = 0, D 2 = h. For h = this 
periodic orbit is replaced by the homoclinic orbit T. When h > 0, the orbits 
with negative D 2 are again periodic or quasiperiodic. On the other hand, the 
orbits with positive D 2 cross into the u 2 > u 2)S region. Their behavior may be 
rather complicated since these may hit the upper boundary of the billiard which 
is slanted. For sufficiently small h, such trajectories closely follow W+ S (P), the 
upper branches of the stable and unstable manifolds of P. Hence, the global 
structure of W+ S (P) determines their behavior. 

We say that W+ a (P) exhibit simple dynamics (SD) if these manifolds in- 
tersect the upper boundary at a monotonically increasing sequence of u until 
exiting the corner region, see Fig. [6] and first column of Fig [7} Otherwise, we 
say that their behavior is complicated, see second and third columns of Fig [7} 
We see that in the complicated cases the manifolds return back to the corner 
region, possibly hitting both boundaries, possibly hitting the boundaries at 

8 Namely for small enough II,— h. 

9 The same behavior extends to all energies in — \{u2, s ) 2 < h < 0, but here we are 
concerned with the near saddle-center behavior. 
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some U\ < U\ s . We describe next the behavior of W™' S (P) and of trajectories 
in their neighborhood near a slightly slanted lower boundary, and discuss the 
behavior of the trajectories at 112 > M2, s in section [6j 

We should note that the behavior of trajectories with a larger oscillatory 
component (e.g. trajectories starting near u 2 , s with I\ > h CTit _^) is expected 
to be complicated as well: almost always such trajectories eventually hit the 
upper ray and we expect that this reflection would, almost always, destroy the 
integrability. Notice that such reflections, induced by the geometry, supply an 
alternative route for the creation of reactions which is unrelated to the local 
structure near P. We will return to this point in the discussion. 



5.2 Non-integrable behavior 

We establish next that for small enough nonzero 9 the motion in the limit 
system is chaotic for a range of energies. We then prove that similar behavior 
occurs for sufficiently small e, c. More precisely, we prove below the existence 
of transverse homoclinic orbits to a saddle periodic orbit (Poincare homoclinic 
orbits). To this aim we use some of the ideas developed in [30] [31] where 
the center-saddle case was analyzed (the first discussion of the problem was 
in [29]). The chaotic nature of the motion follows from this result: it is well 
known that near such homoclinic orbits there is an invariant subset which is 
described on some cross-section by a transitive Markov chain, so in particular, 
this set contains a countable set of saddle periodic orbits, almost periodic 
orbits, etc. [361 EZJ- 

Theorem 5.2 (Complicated Dynamics I) If\9\ is nonzero and sufficiently 
small, then, for e = c = 0, the system (3.3\ ) admits the following propertit 



ties: 



1. The lower branches of the stable and unstable separatrices of P are split; 

2. There is a critical energy value h°, depending on the geometrical pa- 
rameters (9,u s , X,u) , satisfying < h° < h^u-^, such that at energies 
h G (h°, harit-y) the Lyapunov periodic orbit 7/, has two transverse ho- 
moclinic orbits. 

3. At the energy level H = h°(9,u s , X,u) the flow has a tangent homoclinic 
orbit to the related Lyapunov periodic orbit. 
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4- For < h < h° the lower branches of the separatrices do not admit 
sim , pl^\ homoclinic orbits to jh- 



Proof. 1. We construct, as in the proof of Proposition 5.1 the global map 



that is defined by the reflection law from the slightly slanted bottom wall. Let 



\i = tan# and s = (1/a/1 + /i 2 = cos 9, [ij yl + fi 2 = sm9) T be, as above, a 
unit vector defining the lower boundary of the corner. Then, for sufficiently 
small 8, the 3-plane given by u 2 = {iu\ provides a cross-section to the linear 
flow. Recall that the projections of the lower branches of the local stable 
and unstable manifolds of P, W™ ,S (P), are straight lines. Their extensions 
intersect the bottom walj"] at the points: m s = («i )S , 0, X(u 2>s — 

and m u = fiu^ s , 0, — \(u 2 ^ s — /-twi.s))- As before, we consider some small 
cross-sections N s , N u near these points. They are again foliated into 2-disks. 
The coordinates on these cross-sections are (ui,Vi), where here u 2 = [mi, 
and v 2 is expressed from H : v 2 = ±^j2h — 2Ii + A 2 (/iWi — u 2)S ) 2 . The sign 
" + " corresponds to N% and the sign " — " to N%. The restriction of the 
2-form dv\ A du\ + dv 2 A du 2 to (and similarly for N%) is now given as 
dvi A dui + fidv 2 A dv-i, where v 2 is taken from H. Applying the reflection law 



from the slanted wall boundary (4.11) to values (ui,Vi) G iV^, we obtain the 



gluing symplectic map Sh '■ — > (w.r.t. the non-trivial symplectic form, 
see Appendix): 

S h : («!, v t ) i y («!, ( 1 -/^>i + 2 ^2 ) = ^ %i cog 2Q + ^ gin 2Q) (5 3) 

1 + [l l 

This gluing map is defined for trajectories pointing down so that —fivi+v 2 < 0. 
Namely, for sufficiently small we require that v 2 < 0, so the sign in front of 
the radical for v 2 should be negative. 

We now prove the first assertion of the theorem. We show that the reflection 
of the unstable branch (W™(P)) from the lower boundary (the Sg-image of 
the point m u ) does not coincide with the stable branch (W0_(P)) intersection 
with this boundary (the point m s ). For So we should set u\ = Ui )S , v% = 0, 
v 2 = — A(u2, s — ti>U\ yS ). Then we get for /x ^ : 

ui = ui )S , vi = — — - = -2/iAw 2)S (l + 0{fjL)) t 0, 

1 + fi z 



10 Orbits that reflect only once from the corner boundaries. 
11 Recall that since 9 is small, U\ s and u 2 s are positive. 



19 



and therefore the stable and unstable manifolds of P are split - they are of the 
order A|//| apart. 

2-3. Now, let us fix h > 0. We first show that the intersection of Wi(jh) 
(the lower branch of the stable manifold of the Lyapunov periodic orbit 7^) 
with Njl is an ellipse. We then show that the S^-image of the intersection of the 
unstable manifold branch W^( , jh) with N% is another ellipse. We then prove 
that these two ellipses intersect as the energy is increased beyond a critical 
energy. Recall that for < h < h cr it-^ the intersections of Wl^^h) with N s,u 



are bounded away from the cornei 12 and occur, for sufficiently small with 
a vertical velocity which is strictly bounded away from 0. We consider here 
only such h values. 

First, note that the 2-dimensional cylinders W / I'"(7/ t ) are given by solutions 
of the two equations: h(ui, v\) = vf/2 + u 2 (ui — ui tS ) 2 /2 = h, t> 2 = ±A(m 2jS — 
U2), where the +/— signs correspond to the stable /unstable manifolds. The 
intersection of Wt^h) with N s (with coordinates (ui,Vi,v 2 )) is the closed 
curve Ii(ui,Vi) — h, v 2 — A(m 2jS — fJ-Ui), its projection onto the (u±, i>i)-piane 
is an ellipse. Similarly, W u (jh) intersects N u (with coordinates (ui,Vi,v 2 )) 
along the closed curve I\{ui,vi) = h, t> 2 = — A(w 2jS — fiUi), and its projection 
onto the {u\, t>i)-plane is an ellipse as well. 

To find the intersection of the Sh- image of the trace of W^^h) with the 
trace of Wt{lh), h is convenient to use the action-angle (here simply polar) 
coordinates {h,p) on Nfr and (Ii,ip) on N% : V\ = \f2Ly costp, u x = u^s + 



\/2Ii/uj 2 sin p, vi = v 2/i cos ip, u\ = u\ jS + y 2/i/w 2 sin^. Using eq. (5.3) 
the equation Sh{v>i,V\) = (ui,Vi) becomes: 

Ui jS + a/2/i/w 2 sin p = u\ yS + \j2h/uj 2 sinip, 

(1 — fi 2 )\/2h cos ip — 2//A[u2,s — A*(«i iS + \j2h/uj 2 simp 

1 + 



2h cos ip 



The first equation implies that either p = ip or p = 7T — ip. In the first case we 
obtain the following equation for p: 

fxV2h(cos p — — smp) = —X(u 2 s — pM\ s ), 

CO 



12 more precisely, these are bounded away from the 2-plane (ux = 0,U2 = /iitj = 0,^1,^2) 
which corresponds to the corner in the phase space. 
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which has no solutions for sufficiently small In the second case the equation 
for if becomes: 

2h ( cosy? — fi 2 — simp J — — /^A(«2,s — /^l.s)- (5.4) 



Defining 

uj 2 X 2 {u 2 ,s - ^u lyS ) 2 {Xu 2yS ) 2 , N 
a= 2(u* + fi W) =—2- + °M> 

we see that for a fixed /i eq. (5.4) has no solutions for < h < afi 2 , has a unique 



solution at cr/x , and has two solutions for h G (c/x , hcrit-y) ■ Indeed, notice 
that the S^-image of the trace of W^ijh) is given in the («i,i>i) coordinates 
by 

Ul = ^ . A (5.5) 

i>i = cos 28 + -0 2 sin 2$ = — Au 2 s sin 26* + t>i cos 20 + fiXui sin 26 1 , 

that is, it is a linear transformation of the ellipse, so it is also an ellipse. It 
follows from the geometry of intersecting ellipses that for h > a^i 2 these two 
solutions correspond to two transversal homoclinic orbits to the Lyapunov 
periodic orbit 7^, and that at the intermediate case: 

h° = h°(6, u„ A, u) = a/1 2 = ( A ^an0) 2 + Q ^ n q3) ^ 

the unique solution corresponds to a tangent homoclinic orbit of 7^ (it corre- 
sponds to the outer tangency of the ellipses). 

4. Since for h < h° the two ellipses do not intersect, and since in the 
corner interior the flow is integrable, the lower branches do not admit simple 
homoclinic orbits at such energies. In general, it is still possible that conse- 
quent reflections of the extensions of Wt^^h) from the corner boundaries will 
produce homoclinic orbits. ■ 

Next, we assert that similar behavior appears when c and e are sufficiently 
small and the billiard-like potential indeed limits to the impact flow. More 
precisely, in |28j we prove that away from a small boundary layer near the 
billiard boundary, regular reflections of impact flows are C r close to the cor- 
responding reflection-like segments of the smooth Hamiltonians. The smooth 
Hamiltonians are assumed to be in the standard mechanical form with poten- 
tials that are a sum of a smooth potential and a family of billiard-like potentials 
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(satisfying assumptions I- IV of |27J). Additionally, to obtain the correct im- 
pact limit, it is assumed that the value of the full potential along the billiard 
boundary is strictly positive. Then, for positive e the Hill region of the smooth 
flow lies within the billiard region and, near regular impacts, limits to it as 
e — > + . Under these conditions it is proved in [28 J that for sufficiently small 
e the smooth reflection from the Hill region boundary limits to the billiard 
reflection law. Here, to comply with this latter condition, we assume hereafter 
that the potential along the corner rays in At'- 

Rp,L = {{ui,u 2 )\ (u 2 = m tan 9) U (u 2 = «itan(/? + 0)) n (\u h2 \ < L)} (5.7) 

is strictly positive. Namely, we assume that there exist positive constants 
e , B\ such that for all < e < e : 

u 2 A 2 \ 

bV b (R. e u 1 e) + —(u 1 -u 1 , s ) 2 -—(u 2 -u 2 , s ) 2 ) >B X . (5.8) 



For the power law repulsion law (3.5), the first term is infinite, so the inequal- 



ity is satisfied for any b > 0. When the diatomic repulsion is modeled by a 



bounded potential (e.g. (3.4)), b should be taken to be sufficiently large so 
that (5.8) holds. Such an assumption is natural in the chemical-reaction con- 
text: the nuclear potential energy associated with small diatomic distances is 
much larger than the barrier energy. Thus, such an assumption is satisfied 
by adequate models of the PES of triatomic reactions. Then, by [28], we can 
establish: 

Theorem 5.3 (Complicated Dynamics - I - smooth case) Assume that 
Vb(q; e) is a billiard-like potential family limiting to the billiard in the f3-wedge, 



Vf ar fi e id(q) is bounded in the C r topology in the corner region, and (5.8) is 



satisfied. Then, for sufficiently small e,c and 9 max , for all < \9\ < 9 r 



the system (3.3) admits the same properties as listed in Theorem 5.2. There, 
the critical energy levels h°, h crit -^ need to be replaced by a family of (e, c) 
dependent critical energies h e ' c and by h^". Moreover, the critical values 
at which the tangent homoclinic bifurcation occurs depend smoothly on (e, c) 
and approach the limiting bifurcation energy as (e, c) are decreased to zero: 
h £ ' c ->■ h°. 



Proof: Notice that in the proof of Theorem 5.2 for the limit system we 
may replace the cross sections N s with any other locally transverse cross- 
section to W™(P) in the interior of the domain. Consider for example such a 
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local cross section at E" 2 = {(u, v) : u 2 = u 2 < (u 2)S + tan#)/2}. The first 
intersection of the extensions of W^ s {^h) with £ U2 are always transverse for the 
limit system. S" 2 also provides, under specified conditions, a locally transverse 
cross section to the image of the first reflection of the extension of W^i'Jh) from 
the lower boundary. Indeed, for any fixed u 2 one can choose sufficiently small 
\9\ and fix an energy level h U2 (8) so that for all h G [0,h U2 (6)] this section 
is transverse to these orbits. Moreover, h U2 (6) — > /i cr ^_ 7 monotonically as 
\9\ — > 0. Thus, the proof of Theorem 5.2 applies, in the limit system, to the 
traces of W^^h) on E" 2 for all h G [0, h U2 (8)}. Namely, for small \9\, the trace 
which corresponds to the first intersection of the extension of Wt{lh) with S" 2 
and the image of the corresponding trace of W^^h) after its first reflection 
from the lower boundary intersect transversely at h > h°, do not intersect 
at h < h° and are tangent at h — h°. Notice that for sufficiently small \9\, 
h U2 {9) is of order one whereas the ^-dependent homoclinic bifurcation value 



h° is small (see 5.6), namely h° G (0, h U2 (6)). Also, recall that in the proof 



of 5.2 \9\ is assumed to be sufficiently small so that the reflection of W™("fh) 
from the lower boundary is a regular reflection: it occurs with strictly negative 
v 2 . Let 9 max be sufficiently small so that a) Theorem 5.2 applies, b) at 6 max 
the homoclinic bifurcation occurs at an energy at which X" 2 is transverse 
as explained above, and c) this energy is smaller than the nuclear diatomic 
repulsion energy. Namely, at ±9 max , h° < /fe- 7 : = min(/i U2 (±0 max ), Bx). 
This last condition implies that for sufficiently small (c, e), the Hill region of 
(3.3) near the impact point (mi ;S , u\_ s tan#) limits to the lower ray boundary 
of the corner for all h G [0,/fe- 7 ]- Thus, by [SS]. for all \0\ < 9 max and 
h G [0,/i c ^f.J, for sufficiently small (c, e), the stable trace and the image of 
the unstable trace under the smooth flow are C r close to the corresponding 
traces of the impact flow. Hence, the result is established. ■ 

Figure [6] demonstrates that for small \9\ and small e (at c = 0) the lower 
branch of the unstable manifold of P indeed exhibits complicated behavior. 



5.3 Elliptic islands 

Theorem |5.2| implies that for any small fixed fi = tan#, the stable and the un- 
stable manifolds of the Lyapunov periodic orbit 7^ with h° = <Jfi 2 ~ (Am2, s tan 9) 2 /2 
undergo a non-degenerate tangent homoclinic bifurcation. Such a bifurcation 
leads, in particular, to the creation of elliptic periodic orbits near the tangent 
homoclinic orbit: 
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Theorem 5.4 (Complicated Dynamics II) At the semi-interval h > h° , 



for sufficiently small 9, the limit system (3.3) £=c= o has a countable set of 
h— intervals A n accumulating at h = h 0+ such that for h G A n the limit 
flow has generic elliptic periodic orbits with their period tending to infinity as 
n — )■ oo. 

Proof. We show that at e = one may construct a smooth symplectic 
return map to a section close to the Lyapunov periodic orbit 7^0. Then, the 
classical results regarding the emergence of elliptic islands near homoclinic 
tangencies are applicable. 

In order to formulate the classical results precisely, let us recall some details 
[381 139| EOJ HU S2] • Suppose a family of smooth symplectic diffeomorphisms 
fc acts on a symplectic two-dimensional manifold (M,cj). Let p be a saddle 
fixed point for any parameter value c near c = with eigenvalues < A(c) < 
1, A _1 (c). Assume that f 5 undergoes a generic homoclinic bifurcation; at c = 
the stable and unstable manifolds of p have a quadratic intersection at some 
point q (different than p), creating a tangent homoclinic orbit T. Furthermore, 
the unfolding f 5 is generic, namely, in a small neighborhood of T there are no 
simple homoclinic orbits to p at c < and there are two transversal homoclinic 
orbits at c > (or vice versa). Near p one can choose symplectic Darboux 
coordinates (x, y) in which the 2-form u casts as u = dx A dy. The local 
representation of the maps for any sufficiently small c is 

Xi = [A + fi(x,y)]x, yx = [A -1 + gi(x,y)]y, 

with /i,<7i being of the first order at (0,0). Recall that points of V tend to p 
as n — > ±00 (n is the number of the iterations of /). Choose two points q + , g_ 
of T such that q + = (x + ,0) belongs to the segment x > 0,y = (a piece of 
IV s ) and g_ = (0,y_) belongs to the segment of W u , which we assume to be 
y > 0. Then one may take two small neighborhoods V s of q + and V u of g_ 
which belong both to the coordinate chart. Then it is clear that points from 
the semi-neighborhood y > of V s can reach under iterations of / a semi- 
neighborhood x > of V u for all iterations n > n > 1 (n depends on /, on 
the locations of q + , q- and on the size of V s , V u ). The quadratic tangency of W s 
and W u along V means in these coordinates that there is some positive integer 
such that the global map f N : V u — y V s can be written in V u in the form 

x — x + = ax+b{y—yJ) + - ■ ■ ,y = —b^ l x+\Ax 2 +Bx{y—y^) + \C{y—y^) 2j \ 

with nonzero C. It is known that if b > 0, C < then, at c = 0, there are no 
orbits that stay forever in V s U V u under the action of the first return maps 
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/" ° f N (except of T). Suppose, to be definite, that two homoclinic orbits to 
p appear for small c > 0. 



Theorem 5.5 ISffi EH [£TJ In the semi-interval c > there is a 

countable set of c-intervals A n accumulating at c = + such that for c G A n 
the map f E has generic elliptic periodic orbits and their periods tend to infinity 
as n — y oo. 

Recall that the genericity of an elliptic fixed point or a periodic orbit means 
that some of its Birkhoff coefficients in the related normal form do not vanish. 

For our case, the local map f s is obtained as follows. Let (ii,y?) denote 
the symplectic polar coordinates on the (ui,V\) plane (see above). The semi- 
interval ip = 0, < l\ < l\ is a transverse cross-section. We fix h = h° = o~fi 2 , 
and choose, for sufficiently small 6, £, the 3-disk Q : (p — 0, \I\ — I{6)\ < 5, 
\u 2 — Uz, B \ < C> 1^2 1 < Cj where 1(9) = h° is the I\ action of the periodic orbit 
7^0. This Q is a cross-section to orbits that are close to the periodic orbit 7^0. 
The generic family of symplectic maps in the 2-disks Qh is defined by fixing 
the levels H = h in Q, and by defining the bifurcation parameter c = h — h . 
The coordinates on Qh are (^2,^2) and the map in these coordinates is linear 



in «2 — M2,s,^2 : ^ i s obtained from (5.2) if we set t = — in U2(t),V2(t). The 



quarter D 2 < 0, u 2 < u 2yS (recall m 2)S > for the case under consideration) 
corresponds to those orbits of the flow which go from N% to iV^ for positive t. 
The tangent homoclinic orbit T cuts Q^o transversely, when t increases, along 
an infinite sequence of points that lie on the semi-interval D 2 = 0, u 2 < u 2jS , 
v 2 > and accumulate to {u 2yS) 0) - the trace of 7^0 in Q h o. We fix some point 
q + of this sequence. This orbit also cuts transversely the cross-section N^ 
at some point q s + on the trace of W s {^ h o). Since the time of passage between 
these two points of transverse intersections is finite, the flow defines a local 
symplectic map T + : N^ — > Qh°, in some small neighborhoods of q s + and q + 
respectively. 

Similarly, the map T_ is constructed. T_ acts from some neighborhood of 
g_ on the unstable semi-interval D 2 = 0, u 2 < u 2 , s , v 2 < to a neighborhood 
of the trace of T on iV^,- Thus, we have an analytic symplectic map (global 
map) Th = T + o Sh T_ acting from a small neighborhood of g_ to a small 
neighborhood of q + . The global map T^o transforms a small segment of the 
unstable manifold of the saddle fixed point (u 2 ,s, 0) near g_ e Q h o to its image 
in the neighborhood of q + . Since T is a non-degenerate tangent homoclinic 
orbit, this image is quadratically tangent to the stable manifold of the saddle 
fixed point. Moreover, as the mutual position of the two ellipses shows, this 
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segment belongs to the quarter D 2 > 0, v 2 > and u 2 < u 2)S for it. It means 
that we realize the case b > 0, C < of the above mentioned theorem, and we 
get an infinite sequence of intervals in h > h° corresponding to the existence 
of elliptic periodic orbits of the system. ■ 



Theorem 5.6 (Complicated Dynamics II-smooth) Under the same con- 



ditions of Theorem 5.3, for sufficiently small \9\ < 9 max , e and c, at the semi- 



interval h > h E,c {9), the system (3.3) has a countable set of h— intervals A n 
accumulating at h = h e,c (9) + such that for h G A n the limit flow has generic 
elliptic periodic orbits with their period tending to infinity as n -> oo. 

Proof: For sufficiently small \9\, near h°, the gluing map Sh corresponds 



to a regular reflection. Thus, under the same conditions as in Theorem 5.3 
by [2H], for sufficiently small (c, e), the smooth version of the return map to 
Qh, Tfr C , is C r close to Th for all \9\ < 9, max and h G [0,/i^£!:J. Hence, by 



theorem 5.3, for h > h £ ' c (9), where h £ ' c {6) denotes the energy of the homoclinic 



tangent bifurcation of the smooth system, theorem |5.5 may be applied to the 



return map T^' c . Moreover, this return map depends smoothly on (e, c), hence 
the theorem follows. Notice that the classical results are clearly applicable for 
finite e values for which the homoclinic tangency persists. Here, we see that 
one may change the order of the limits, namely, even for arbitrarily small e 
the stability islands appear. ■ 



6 Criteria for simple dynamics 

The behavior near the lower branches of the stable and unstable manifolds of 
P was shown to be integrable when 9 = e = c = and non-integrable for small 
non-zero 9. Indeed, we proved that the stable and unstable manifolds of 7^ 
intersect for a range of energy values, h G (h° m ( AM2 .^ ap61 ) ; h cr i t - 1 ). Similar 
results apply for the upper branches when (/3 + 9) is small. 

Proving analogous results regarding integrability or non-integrability of the 
dynamics near P for general (/?, u s , L, 9, u, A) is a difficult problem. Rather 
then seeking the general solution, we identify cases in which the nature of the 
dynamics can be roughly predicted. More precisely, we classify the behavior 
of the upper/lower branches of the manifolds as follows: 
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SD (Simple dynamics) For small h, the upper (respectively lower) branches 
of the unstable and stable manifolds are reflected to the regions of "no- 
return" monotonically (with a monotone sequence of reflecting points 
exiting the corner region) and thus, near these branches there are no 
recurrent trajectories in the neighborhood of P. 

PCBD (Possibly complicated bounded dynamics) For small h values 
the upper (resp. lower) branches of the manifolds are trapped in the 
upper (resp. lower) corner region: these manifolds cannot exit the cor- 
ner region from the product (resp. reactant) channel. 

CD(Chaotic dynamics) The branches of the manifolds intersect each other 
for a range of energy levels that are close to the barrier energy. 

The first and second rough criteria check whether the manifolds reflect in 
a simple, monotone way out of the corner region or whether they are trapped. 
Clearly simple dynamics imply that no chaotic behavior is possible for the 
corresponding branches. Folding back of the manifolds is a necessary criterion 



for the creation of homoclinic tangles yet it is not sufficient, as the special 



case e = 9 = shows. Hence we call this behavior possibly complicated. 

Notice that when both the upper and lower branches exhibit SD the linear 
structure near the saddle point governs the motion. On the other hand, when 
one set has CD and the other SD we have the "open chaos" scenario in the 
reaction region. In particular, then the reactant/product region is strongly 
asymmetric. When both branches intersect we have the classical double loop 
homoclinic tangles. If additionally the PCBD conditions are satisfied, the 
chaotic motion associated with this homoclinic tangle is limited to the reaction 
region (and then the implications regarding scattering have yet to be explored). 

Next, we show that there are some regions in the parameter space where we 
are able to determine that SD occurs and others where PCBD occurs. Figure 
[8] summarizes some of these results graphically. 

First, we notice that independent of the saddle eigenvalues (i.e. of co,X), 
when the upper and/or lower branches of W± S (P) do not reflect from the 
upper (respectively lower) boundary ray we have simple dynamics: 

Theorem 6.1 (Simple Dynamics I) For sufficiently small e,c, if (3 + 9 > 
it/2 or ifu\ >s tan(/3+#) > L, there exists h** (ui )S , (3+9, L\ e, c) > 0, such that for 



13 Trapping of the manifolds implies the folding back of them but the opposite implication 
is not true in general - the lobes of homoclinic tangles may extend to infinity. 
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h G [0, h**), H / "' s (7/j) ; the upper branches of the stable and unstable manifolds 
of the Lyapunov orbit 7^ exit the corner region without intersecting each other. 
If tti jS tan(0) < — L, fork < h**(ui tS ,—9,L;e,c) the same statement applies to 
the lower branches, W^ s {^h)- 

Proof: Consider the e = c = case. Here the projection onto the configuration 
space of the extensions of W+ S (P) is a straight vertical segment that intersects 
the upper boundary of Al at (ui,u 2 ) = (ui >s ,L). The corresponding phase 
space intersection points of W+ S (P) with the three dimensional cross section 
Q = {(u, v)\u2 = L} are (u±, Vi, u 2 , v 2 ' s ) = (ui }S , 0, L, ±A|L — U2 )S \), so, these 
points are well separated in phase space. Similarly, if h is not too large (see 
below), the traces of the stable and unstable manifolds of 7/1 on Q are two 
planar ellipses, parallel to the (ui, v 1) plane and separated by the finite distance 
2X\L — u 2>s \ along the v 2 axis, so they do not intersect (this is simply the 
linear behavior). The restriction on h appears from the requirement that the 
manifolds should not hit the upper ray of the corner for any u 2 ^ L. An 
explicit bou nd on h may be thus easily found from setting I\ = h and v 1 = 



in eq. (4.7). If (3 + 9 > ir/2 then h < h* Q l t = \(u liS ) 2 (so that 7 h does not hit 



the upper ray). If U\ jS tan(/3 + 6) > L, we require Ui(t) < Lj tan(/3 + 6) for all 
t, hence h < h** cute = ^- ( tan( -^ +e - ) — Mi, s ) 2 - The same arguments apply to the 
lower branches when 9 < (replacing L by — L and ((3 + 6) by 9). These results 
are concerned with robust properties of trajectories within the corner region 
(with no impacts), hence, by the smooth dependence of trajectories within the 
corner region on parameters these are clearly true for sufficiently small e, c. ■ 

The second geometrical observation is that the projection of W+^i^h) onto 
the configuration space lies within the Hill region's of the energy level h. In 
particular, for P = 70, this Hill region boundaries are the potential zero- 
level lines V a (u) = together with the corresponding boundaries of the corner 
region Al (see below). We then notice that for sufficiently large j these zero 
potential level lines always intersect the corners' boundary at U\ < L, namely 
that W+ S (^h) must fold back inside Al- We thus find the critical value of j 
above which the folding occurs: 

Theorem 6.2 (Possibly Complicated Bounded Dynamics I) For any given 
geometrical parameters (u s , (3 + 9,L) with L > tan(/3 + 6)u^ s and (3 + 9 < n/2, 
there exists <5+ > tan(# + (3) as specified below, such that for all y > 5 + , for 
sufficiently small h and (e,c), the upper branches of the stable and unstable 
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manifolds of the Lyapunov periodic orbit 7^ (W+^^h)) do not exit the corner 
region through the product channel (i.e. with U2 > u 2iS ). 

Proof: Consider the region bounded by the corner and the upper segments of 
the zero level lines of the quadratic potential: 

C + = At H {('«!, W 2 )|(«2 - U 8>2 ) > j\Ux -W s ,l|} ■ (6.1) 

C + is the Hill region for the energy level h = for orbits of the limit system 
belonging to W^ s (jo = P); since the kinetic energy must be non-negative, 
orbits belonging to the energy surface h = must satisfy for all time V(u(t)) < 
0. Hence, at e = c = 0, W+ S (P) is indeed trapped in C + . Similarly, orbits 
belonging to W^' s ( , ~fh) for h > are restricted to reside in the configuration 
space region C\ at which V(u(t)) < h, and initially (eventually for the stable 
manifold) also satisfy u% > u 2>s . 

The shape of the region C + depends on the geometrical parameters. When 
tan(0 + (3) < j, the upper ray of the corner u 2 = U\ tan(/3 + 9) intersects the 
upper zero-potential energy level line {(wi, 1*2)1^2 — u 2>s = j( u i ~ u i,s), u 2 > 
U2, s } at the vertex point (ui,u 2 ) with 

+ U2 ' S " fUl ' S u x , s + ta y + ^~r 2 '% ut = ui tan(/3 + 9). 



1 tan(/3 + 0)-f i,s f-tan(/3 + 0) 
Let 5c = t be the eigenvalue ratio for which L = max(tt|, u 2 



5+(L,u s ,ttm((3 + 9)) 



tan(/3 + 9) + f or (3 + 9<ti/A 



tan(/3 + 0) + y^g; for /3 + ^>-/4 

(6-2) 

Notice that for ^ > 5+ the vertex is inside the corner region, (max(w+, u 2 ) < 
L). Fig. [I] presents the typical dependence of 5+ on 9 and L (blue curves). 

Similarly, for sufficiently small h>0, when ^ > 8£, the right boundary of 
the region (the upper part of the h Hill region), intersects the upper corner 
boundary at some finite value {u\ (h) , u 2 (h) = w*(/i)tan(/3 + 9)). Then, the 
shape of the upper part of the Hill region is triangular like, and orbits in its 
upper part have U\{t) < u\{h) < L. Thus, for small h, trajectories belonging to 
W+ S (jh) cannot escape the corner region with u 2 > u 2jS , namely, generically, 
the behavior is not simple. Finally, for sufficiently small (e, c), the restricted 
Hill regions are deformed into Hill regions with the same basic property: when 
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t > for sufficiently small h ^ 0, for w 2 > w 2)S , the Mi coordinate is limited 
by some u\{t) < u\{h,e,c) < L. ■ 

Note that the assumption that c is small is equivalent to requiring that 
the far-field terms remain small in Al - otherwise, these far-field terms may 
indeed change the shape of the Hill regions in Al- 

A similar statement for the lower branches is: 

Theorem 6.3 (Possibly Complicated Bounded Dynamics II) For any 

given geometrical parameters (u s ,8 < 0,L) there exists 5~ > — tan# as spec- 
ified below, such that for all j > 5~ , for sufficiently small h and (e,c), the 
lower branches of the stable and unstable manifolds of the Lyapunov periodic 
orbit 7^ (W"' s (7/ i ) y ) do not exit the corner region through the reactant channel 
(i.e. with u 2 < u 2>s ). 



Proof: As in Theorem 6J2, we find the lower part of the Hill region for h = 
£ = c = 0: 

C~ = An {(ux,u 2 )\(u 2 - U s>2 ) < °^\ui - U Sj i\j . (6.3) 

If tan# > —j, the lower ray of the corner u 2 = uit&nO intersects the lower- 
right zero-potential energy level line {(ui,u 2 )\u 2 — u 2<s = —j(ui — u\ iS ),u 2 < 
u 2,s} at the point {u^.u^ ) with 

u 2>s + fwi, s , u 2)S - u x>s tan 9 _ 

m, = — — = Mi s H — - — , Un = Ui tan v. 

1 tan 9 + f * tan 9 + f ' 2 1 

Setting, for 9 < 0, 

f -tan 9+ ^- tane "M for _# <7r / 4 
i ~tang+ _ L/tan0 _ u ' a for -6 > tt/4 

we obtain that for all j > 5~(L, u s , tan 6) the lower branches are trapped in 
the corner region, see Fig. [8] for the typical dependence of 5~ on 9 and L (red 
curves). Repeating the same arguments as for the upper branches (Theorem 



6.2) the theorem is established. ■ 
Notice that for 9 > 0, for all |, u{ < namely the Hill region is always 
bounded. Moreover, the manifolds in this case are reflected towards the corner 
region; Recall that the first hit of the lower boundary of trajectories belonging 



to W™(P) is with v 1 = 0,v 2 < 0, hence, by (4.11), vf < 0. We conclude 
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that these trajectories reflect towards the corner region and hence the flow is 
not simple for e — c — h — 0. Using continuity and the regular e limit at 
regular reflection we see that a similar statement applies to the small h, e, c 
case. Hence, we conclude: 

Corollary 6.4 (Possibly Complicated Bounded Dynamics III) When 9 > 
the lower branches of the manifolds fold back and the flow is possibly com- 
plicated. 

Notice that when e > 0, the appearance of closed level sets of the potential 
function (as established in the above theorems and corollary), implies that 
these regions also contain a minimum point of the smooth potential. Thus, we 
see that for large co/X, even though the limit system has a single saddle fixed 
point in the corner, the smooth system has several fixed points and some of 
them are stable. 

Finally, Figs. 7|9fl~0 suggest that for 9 < and sufficiently large A simple 



dynamics are realized: 

Conjecture 6.5 (Simple dynamics II) For any fixed geometrical parame- 
ter (u s ,9, (3, L) satisfying —f3 < 9 < 0, for a fixed u and sufficiently large A 
(small S — j), for small h and sufficiently small {e,c), both the upper and 
lower branches of the stable and unstable manifolds of the Lyapunov orbit 7^ 
C^± S {lh) ) have simple behavior - these exit the corner region through the prod- 
uct (respectively reactant) channel without intersecting each other in the corner 
domain. 



Supporting evidence: One may start by proving the claim for h = e = c = 0, 
proving that trajectories belonging to W+ S (P) reflect only a finite number of 
regular reflections before exiting the corner region through the right side of 
a box of size L. Then, the claim follows by the smooth dependence of the 
manifolds on parameters (for sufficiently small h, c) and by the closeness of 
the limit system to the smooth system at regular reflections (for sufficiently 
small e). Let (u(t),v(t)) G W™{P) hit the upper corner ray at times t iy i > 1, 
so that t± is its first hit. Numerical simulations show that as A is increased the 
sequence of Mi(tj) is indeed monotonically increasing. Moreover, these suggest 
that the gaps M 1 (t i+1 ) — Wi(tj) approach a constant value: in between hits the 
v 1 velocity remains essentially constant since the time between reflections is 
roughly of order 1/A 2 whereas the changes in v 1 are on the much longer time 
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scale 1/u 2 . Proving the above statements requires quite elaborate calculations 
of the asymptotic behavior at small 5, calculation that go beyond the scope of 
this manuscript. 

Figures [7] and [9] demonstrate that the three different behaviors may be 
realized in the singular limit (e = in Figure [7]) and in a smooth case (e = 0.01 
in Figures [9]) when we change A and keep all other parameters fixed; The left 
columns show simple dynamics in both the upper and lower branches, the 
middle ones simple dynamics in the upper branch and possibly complicated 
dynamics in the lower branch and the right panels show possibly complicated 
dynamics in both branches. Fig [8] shows that these findings are consistent with 



the 8* bounds for j. Figure 10 shows the regularization effect that is achieved 



when e is increased. Notice that even without introducing far-field potential 
terms, the geometrical potential level curves are reminiscent of Fig[l]D and are 
similar to other PES appearing in the Chemistry literature PQ-[9]. 



7 Discussion 

The stable and unstable manifolds of unstable periodic orbits with energies 
that are slightly above the barriers' energy divide the initial iso-energetic 
phase space region of incoming trajectories to reacting vs non-reacting regions 
[3 El EE2] • The structure of the manifolds is called simple if these manifolds 
do not intersect each other and simply extend to the reactant and product 
channels. Then, the phase-space transition state theory provides an accurate 
description of the transition rates from reactants to products. On the other 
hand, if the manifolds intersect each other or fold back into the reaction region, 
there is no cross-section which is crossed only once by all incoming trajecto- 
ries. Then, the main assumption underlying the transition state theory fails 

[3 Sua EE] . 

By introducing a geometrical model for the reaction dynamics we find con- 
ditions under which the manifolds structure is simple and conditions under 
which it is complicated. Three qualitative observations emerge. First, we 
proved that a homoclinic bifurcation occurs when the manifolds are close to 
being perpendicular to one of the corner rays. Then, there are intervals of 
energies at which stable periodic triatomic configurations emerge. Second, we 
found that when the projection of the unstable eigenspace to the configura- 
tion space intersects the lower corner ray in an obtuse angle (so 9 < 0), and, 
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additionally, the saddle-center expansion rate is much larger than the oscilla- 



tion frequency, the manifolds geometry is simple (Conjecture 6.5). Third, we 
established that provided the unstable eigenspace direction intersects the cor- 
ner rays, the manifolds are trapped for sufficiently large oscillation frequency 
(Theorems 6.2|6.3 ). 



We expect that similar qualitative statements may be formulated when non- 
linear normal form and higher dimensional extensions are included. Namely, 
that conditions under which phase-space transition state theory is adequate for 
describing the reaction for energies that are close to the saddle energy may be 
found by a similar methodology. More generally, while the effects of non-linear 
terms in the reaction region have been suppressed here for simplicity of pre- 
sentation (by taking only quadratic terms of the integrable normal form and 
by considering small c), we trust that the main principles that were discovered 
hold for the nonlinear case as well; Indeed, once the fixed point P has the 
saddle-center structure, its stable and unstable manifolds may be computed. 
With a general non-linear smooth potential their projection to the configura- 
tion space will appear as curved lines that may or may not intersect the corner 
region boundary. The analysis presented here applies to the case where, in the 
limit system, the manifolds do hit the boundary and reflect back. Then, we 
expect to find similar behavior in terms of the Hill region dependence on the 
ratio between the oscillatory and the hyperbolic eigenvalues. 

Finding the implications of the above observations for specific chemical 
reactions is an interesting and challenging endeavor - we believe the tools 
developed here may shed some light regarding the governing parameters. 

Finally, we note that the analysis presented here applies to the traditional 
energy regime in which the behavior near the barrier is examined. The behav- 
ior for larger energies, or equivalently, for reactions that do not have a barrier, 
is expected to be quite different and will be described elsewhere. Indeed, we 
hold that there are two distinct mechanisms that give rise to the observed 
sensitive dependence of the reaction rates on the energy and the initial con- 
ditions: those associated with the complicated structure of the manifolds as 
discussed here and those associated with the corner geometry of the nearly- 
billiard Hamiltonian. 
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Appendix: The gluing map is symplectic. 

Let us check, for the reader convenience, that the map Sh '■ N% — > defined 
by the reflection law is a symplectic map w.r.t. restrictions of the main 2-form 
on and N%, respectively. We shall verify it for the lower wall supposing 
H = tan# is finite. We work in small neighborhoods of the points m„ and 
m s being the intersection points for the lower branches of unstable and stable 
manifolds of the equilibrium P. These cross-sections belong to the 3-plane 
given by the relation ui = fiu±, the restrictions of 2-form dv\ A du\ + dv2 A du2 
to N%, Njl are the following: 

Cj = dvi A dui + tan^rf-0 2 A du\ = (1 — tan 6^)dv\ A dui, 
u = dvi A du\ + tan#<it>2 A du\ = (1 — tan6^)<ifi A du\. 

The symplecticity condition for Sh means, as is known S^u = u [13]. Using 
the relations Ui = u\, V\ = , 0icos2# + -0 2 sin2^ and expression for v 2 on iV^ 4 
one gets: 

Stuj = f 1 - tang ! 1 cos2 ^ + ^ 2 sm26 j [dv lC os26 + dv 2 sm29] A dui = 
h V t)isin2g-{) 2 cos2gy L J 



v\ cos 29 + v 2 sin 29 

1 — tan 9- 

v\ sm 29 — v 2 cos 29 

1 — tan 9^- ) diii A du\ = Cj. 



cos29 - sin 29 — 

V2 



dvi A dii\ 



v 2 . 

Thus we get that the Poincare map o Sh is also symplectic. 
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(a) 



(b) 



Figure 1: Contour lines of the effective potential for the Hi + H reaction 
in the (a) relative positions coordinates V(ri,r 2 ) (b) mass scaled coordinates 
V r (qi, qi). The allowed region of motion lies within the dense contour lines (red 
lines) that correspond to the strong diatomic repulsion at small distances. In 
the (r±, r 2 ) plane this region is essentially the positive quadrant whereas in the 
(q±, q 2 ) plane it is confined to a (3- wedge/corner: a two-dimensional wedge with 
a corner angle = 60°. The appearance of a single saddle point in the corner 
region is apparent. The plotted potential is of the LEPS form, see [H|2j|3] and 
references therein. Notably, the /3-wedge feature always appears in the mass 
scaled coordinates, independent of the exact form of the potential. 
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Figure 2: The saddle in a corner geometry in the configuration space. The 
potential level sets near the saddle are schematically drawn; The zero level 
lines are the thick black lines whereas the thin black curves are the nearly- 
zero potential level lines. The projection of the center eigen space is denoted 
by the red line and that of the stable and unstable eigen spaces is denoted by 
the green line. The top (bottom) panels show the geometry when 9 is positive 
(negative). The left (right) panels show the geometry in the mass scaled (qi, #2) 
(the rotated (ui, U2)) coordinates. The region Al is indicated by grey shading. 
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Figure 3: The phase space structure of the linear system. Shaded region with 
dashed level curves correspond to non-reacting orbits having D 2 < (these 
trajectories cannot pass the dividing configuration-space section u 2 = U2, s )- 
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Figure 4: The return map for 6 = 0. The green triangle on N s is mapped 
by the smooth flow to the green triangle on N u . Its reflection (the image) 
is the black triangle belonging to N s . Here, it belongs to the same level sets 
Ii,D 2 yet its velocity v 2 changes sign. 
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Figure 5: The return map for small 9. As in Fig 4, the image of the green 
triangle of N s under the flow and the gluing map is the black triangle. Here, 
the velocities (fi,^) and the corresponding level sets Ii,D 2 are changed due 
to the reflection. 
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(a) 



(b) 





(c) 



(d) 



Figure 6: Projections of the unstable manifold of P onto the configuration 
space at small 9 values, (a) 6 = 0.1 (b) 6 = 0.01 (c) 9 = —0.1 (d) Zoom-in 
of (b). Thick blue (red) line: W$(P) (W™(P)). Thin colored lines are the 
level sets of the geometric potential, with the exponential repulsion form (Eq. 



(3.4)). Here, b = 10, e = 0.01, u s = (2.5, 1.5), = n/3,u 2 = 1, A 2 = 3. 
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Figure 7: Simple and complicated dynamics for the limit system (e = c = 0). 
The projections of W+(P) (blue) and W^(P) (red) to the configuration space 
(top), (ui,Vi) space (middle) and (u2,V2) space (bottom) are presented for 
three different A values: left, middle a_nd right columns have A 2 = 2., 1.4, 0.5 
respectively. In all figures u 2 = l,u s = (2.5, 0.3), (5 = n/3,9 = —n/6 and 
L = 6. The colored lines show the level curves of the quadratic potential 
V a [u) (top), the center constant h(ui, vi) (middle) and the hyperbolic constant 
D2{u2,V2) (bottom). The magenta rays on the top panel indicate the impact 
surfaces and the magenta squares indicated the location of the stagnation point 



Figure 8: The dependence of the critical values of w/A on ^ and L. The 
critical values 5f (see Theorems 6.2|6.3 ) for the upper (blue) and lower (red) 
branches of the manifolds are plotted for increasing L values. When u/X > 5 ± 
the upper/lower branches of P are trapped in the corner region. Here, 
(3 = 7r/3,L = 6, 26, 46, .., 206, u s = (2.5,0.3). The convergence at large L 
values is apparent. Simple behavior is expected to appear at small u>/\, below 
both critical curves. The values of tu/\ for the three A values of Fig. [7] are indi- 
cated by squares (lowest, middle, upper squares correspond to left, middle, right 
columns of Fig [7]) . 
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Figure 9: Simple and complicated dynamics for the smooth system. The 
same parameters, initial conditions and projections of Fig. [7] are used, yet, 
the impact flow is replaced by a smooth exponential potential along the rays 
(Eqs. ( 3.3|3.4 ) with b — 10, £ — 0.01^ The colored lines on the top are the 
level curves of the smooth potential V a (q) + bVt(q, e). 



Figure 10: Simple dynamics for the smooth system: the effect of smoothing. 
The same parameters, initial conditions and projection of Fig. 6 top left 
panel are used with the smoothing parameter increased to £ — 0.05 (left) and 
e = 0.1 (right). The smooth potential level curves still follow those of the 
linear potential quite closely in the region explored by W±(P) . 
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